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Abstract 

We propose a formalism inspired by matrix models to compute open and closed topo- 
logical string amplitudes in the B-model on toric Calabi-Yau manifolds. We find closed 
expressions for various open string amplitudes beyond the disk, and in particular we write 
down the annulus amplitude in terms of theta functions on a Riemann surface. We test 
these ideas on local curves and local surfaces, providing in this way generating functionals 
for open Gromov-Witten invariants in the spirit of mirror symmetry. In the case of local 
curves, we study the open string sector near the critical point which leads to 2d gravity, 
and we show that toric D-branes become FZZT branes in a double-scaling limit. We 
use this connection to compute non-perturbative instanton effects due to D-branes that 
control the large order behavior of topological string theory on these backgrounds. 
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1 Introduction 

Topological strings constitute an important subsector of string theory with various physi- 
cal and mathematical applications, and they have been extensively investigated since they 
were first formulated. This has led to many different ways of computing their amplitudes 
based very often on string dualities. Topological strings come in two types, the A model 
and the B model, which are related by mirror symmetry. The A model provides a physical 
formulation of Gromov-Witten theory, and it can be solved on toric backgrounds by a 
large N duality with Chern-Simons theory. The B model is deeply related to the theory of 
deformation of complex structures, and the amplitudes in the closed string sector satisfy 
holomorphic anomaly equations [H] which can be used to solve the theory in many cases. 

The open string sector of topological string theory is not as well understood as the 
closed one. For example, the open string analogue of the holomorphic anomaly equations 
has not been formulated yet. Most of the explicit results in the open B model refer 
to disk amplitudes, which have been studied in detail in the local Calabi-Yau (CY) case 
PS El EE, EDI El] and more recently in the compact case [70], [65]. In [2] a general framework 
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was proposed to understand the B-model on local CY geometries, but in practice its use 
has been limited to the simplest local model i.e. the topological vertex. 

On top of these practical problems concerning the computation of the amplitudes, 
there are more conceptual problems related to the nature of the string perturbation series. 
Like other string theories, topological strings are defined in a perturbative, genus by genus 
expansion, and we do not have a general definition beyond perturbation theory. Some 
proposals have emerged based on holographic descriptions related to black hole physics 
[63] . and in some simple cases these descriptions have been worked out in detail [6], but 
much remains to be done in order to understand these fundamental issues. 

In this paper we discuss the open string sector of the B model and the nature of the 
topological string perturbation series by relying on ideas and techniques of matrix models. 
Dijkgraaf and Vafa have shown that, in some local Calabi-Yau backgrounds, the B model 
reduces to a matrix model [28], but the backgrounds they considered do not include 
the toric backgrounds with an interesting enumerative geometry. Some of these, like the 
resolved conifold and the An fibrations over P 1 that engineer M = 2 gauge theories, admit 
a matrix model description closely related to Chern-Simons theory [531 E] • However, the 
matrix model that one finds in this way is difficult to solve and its applications have been 
limited. 

A way to remedy this situation is to notice that the 1/N expansion of a matrix model 
can be substantially generalized and reformulated, leading to a geometric construction 
that starts with a Riemann surface and defines recursively a set of open and closed am- 
plitudes on it. This construction emerges in a natural way when one solves the loop 
equations of a matrix model. It was started in [12] and has culminated in the work of 
Eynard and collaborators (see [221 EDI El], and specially [32], where this point of view 
is stated most clearly). The general construction includes as a particular case the 1/N 
expansion of the standard one-matrix model with a polynomial potential. In that case, 
the spectral curve encoding the planar solution can be described by an equation of the 
form, 

n 

y 2 (x) = M 2 (x)a(x), a(x) = JJ(x - Xi ), (1.1) 

t=i 

where M(x) is a polynomial in x. The corresponding closed and open amplitudes ob- 
tained from the general geometric construction of [291 E2] are then the free energies and 
correlation functions, respectively, of the matrix model underlying (11.11) . 

For an arbitrary Riemann surface we should expect that the quantities obtained in 
this formalism are general amplitudes for a chiral boson living on the curve, since the 1/N 
expansion of a matrix model can be reformulated in this way [37| . On the other hand, it 
was shown in [2] that the B model on a general local Calabi-Yau geometry is described as 
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well by a chiral boson living on a Riemann surface. It is then natural to conjecture that 
the geometric amplitudes defined in [29], [2H [32] describe open and closed string amplitudes 
for the local B model, once the appropriate spectral curve has been found. We will see 
in this paper that this is the case also for toric geometries. The spectral curve can be 
obtained from the disk amplitudes, and it is of the form (11 .11) but with a nonpolynomial 
M(x). Using this we give closed B model formulae for various open string amplitudes on 
toric CY threefolds beyond the disk level, like annulus amplitudes and three-holed sphere 
amplitudes, which we test in various examples. In these tests we focus on so-called outer 
brane amplitudes of toric geometries. We expect that our formalism can accommodate 
more general cases, and we offer a preliminary discussion of this issue in the concluding 
section. 

Another advantage of matrix models is that one has a good control of instanton ef- 
fects in the form of eigenvalue tunneling [22| [68] . This has been used to understand the 
properties of the string perturbation series in the double-scaled models which describe 
non-critical strings, by using the standard connection between instanton effects and large 
order behavior in perturbation theory [49] . In [18] it was found that topological strings 
on local curves provide a realization of 2d gravity near a critical point at small radius 
in Kahler moduli space. This was done through a detailed analysis of the closed string 
sector. Here we analyze the open string sector of these models to test our general matrix 
model formalism, and we show that, in the double scaling limit, toric D-branes become 
the FZZT branes of 2d gravity. This allows us to identify computable, non-perturbative 
instanton effects due to D-branes which are responsible for the large order behavior of 
topological string perturbation theory away from criticality. Our analysis suggests that 
in some circumstances one can find Borel-summable string perturbation series. 

The organization of this paper is the following. In section 2 we review open string 
amplitudes in the local setting and the relevant matrix model results, and we present our 
proposal for the computation of open (and closed) string amplitudes in the local B model. 
In section 3 we analyze the example of local curves. We test our proposal and we analyze 
the instanton effects controlling the large order behavior. Section 4 presents further tests 
for two local surfaces: local P 1 x P 1 and local P 2 . In section 5 we conclude with various 
problems and comment on possible developments. 

2 Open string amplitudes and matrix models 

In this section we first review various aspects of topological string theory relevant to this 
paper and then we make our proposal. General recent references on topological strings 
include [541 l69l [59], to which we refer the reader for background material. 
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2.1 General aspects and open Gromov— Witten interpretation 

Let us consider topological string theory on a Calabi-Yau manifold X. To define the 
open string sector in the A model we have to speciy a Lagrangian submanifold £ of X 
|71j . The type-A open topological string theory describes holomorphic maps from open 
Riemann surfaces of genus g and with h holes to the Calabi-Yau X, with Dirichlet 
boundary conditions specified by C. These holomorphic maps are usually called open 
string instantons. The topological sector of an open string instanton is given by two 
different kinds of data: the boundary part and the bulk part. For the bulk part, the 
topological sector is labelled by relative homology classes, since we are requiring the 
boundaries of /*[£ ff /J to end on C Therefore, we will set 

MZ a , h ]=(3eH 2 (X,£)- (2.1) 

To specify the topological sector of the boundary, we will assume that = 1, so that 

H\(C) is generated by a non-trivial one-cycle 7. We then have 

f,[Ci]=wa, Wi G Z, i = l,---,h, (2.2) 

i.e. Wi is the winding number associated to the map / restricted to Cj. We will collect 
these integers into a single /i-uple denoted by w = (wi, • • • , w^). 

The free energy of type-A open topological string theory at fixed genus and boundary 
data w, which we denote by F w g (t), can be computed as a sum over open string instantons 
labelled by the bulk classes: 

FgA t )=J2 F 9,^Q P - ( 2 - 3 ) 

P 

In this equation, the sum is over relative homology classes (3 G H 2 (X, C). The quantities 
Fw,g,p are open Gromov-Witten invariants, and they 'count' in an appropriate sense the 
number of holomorphically embedded Riemann surfaces of genus g in X with Lagrangian 
boundary conditions specified by £, and in the class represented by fl,w. Finally, Q 13 
denotes exp(— J2i PiU), where the index i = 1, • • • , 6 2 {X, C) runs over a basis of H 2 (X, C). 
It is useful to have a generating functional for open string amplitudes at fixed genus 

F 9 (z 1 ,---,^) = ^F^(t)^ 1 "--C- (2-4) 

Wi 

Particularly important cases of these amplitudes are the disk and the annulus amplitudes, 
which correspond to g = and h = 1,2, respectively. We will denote them by W(z) 
and A(zi,z 2 ). These amplitudes have the following spacetime interpretation (THE?]. If 
one considers N D6 branes wrapping the Lagrangian submanifold C, the resulting four- 
dimensional worldbrane theory is an U(N) gauge theory with M = 1 supersymmetry 
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and chiral fields associated to the closed and open moduli t, z. The superpotential for 
these moduli is precisely the disk instanton amplitude W(z), while A(z, z) gives the gauge 
kinetic term. 

It should be possible to compute these open string amplitudes directly in the B-model, 
and then use mirror symmetry to relate them to the amplitudes in the A model, as it 
was done in [3 E] for the disk. A special case where we know the B model answer is 
the situation considered by Dijkgraaf and Vafa in [28]. In the local CY backgrounds of 
[28] the B-model reduces to a matrix model for a Hermitian matrix M, and open string 
amplitudes have the following interpretation. Define the connected /i-point functions, 

W{p x , ■ ■ ■ ,p h ) = <Tr -J— ■ ■ ■ Tr^— > {c) , (2.5) 

which are generating functionals for matrix model correlators 

W{ Pl ,-- - ,p h ) = ^(TrM fcl ■■■TrM kh ) (c) p- kl - 1 ■■■p^ h_1 . (2.6) 

ki>i 

These quantities have a genus expansion of the form, 

oo 

W( Pl , ■ ■ ■ ,p h ) = Y^9 2 r 2+h W 9 ( Pl , ■ ■ ■ ,p h ). (2.7) 

9=0 

If we set Zi = pi i the relation with (12.41) is simply 

W g { Zl ,--- ,z h ) = z 2 1 d Zx ---z 2 h d Zh F g {z 1 ,--- ,z h ), h>2. (2.8) 

The case of g = 0, h = 1 is slightly different. The generating functional is the resolvent 
of the matrix model, and it is related to the disk amplitude as 

ujq(z) = z + z 2 d z W(z). (2.9) 

In the case of a one-matrix model with a polynomial potential V(M), the resolvent is 
always of the form (see for example [251 155] ) 

Mp) = \(V(p) + y(p)), (2.10) 

where y(p) is the spectral curve describing the planar model and can be written as in 
( 11. II) . The points Xi are the branch points of the curve. 

Consider now the case in which A is a toric CY threefold. The mirror CY can be 
written in terms of an algebraic equation, 

F{e u ,e v ) = 0, (2.11) 
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which defines implicitly a function v(u). It was shown in [7] that in these geometries there 
are Lagrangian submanifolds C with the topology of S 1 x C whose disk amplitudes can be 
easily computed in the mirror geometry (12. lip . As shown in [7J [S], they are simply given 
by 



W(u) 



du'v{u'). (2.12) 
There are various subtleties in 



Here, one identifies the open string modulus as z = e 
this procedure. First of all, there are various types of Lagrangian submanifolds in these 
geometries depending on their location in the toric diagram. The choice of parametrization 
of the curve (12. lip depends on this location. Second, the variables appearing in (12.111) 
are not flat coordinates. There is an open/closed mirror map of the form 



U = U(u,t), 



(2.13) 



where the flat open coordinate depends on the closed string moduli t. One has to reexpress 
the superpotential in terms of the U variable in order to extract geometric information. We 
will discuss these issues in more detail in section 4, when we consider concrete examples. 
In any case, the above results suggest to define a resolvent for these local models as 



^o(p) = -(1 + v(p)), 
P 



(2.14) 



The function v(p) has in general the structure 

a(p) + y/a{p) 



v(p) = log 

If we now use the identity, 

log 



c(p) 



(2.15) 



i=i 



a + y/a 



1 , a 2 — a 
- log h tanh 



-i 



(2.16) 



we see, by comparing ( 12.101) and (12.15p . that the role of the spectral curve is played in 
these geometries by the equation 



i \ — i 

y(p) = - tanh 

V 



vn?=i(p- 

a(p) 



Xi 



(2.17) 



This can be written in the form ( ll.ip but with a non-polynomial function M(p): 

2 



M(p) 



tanh 1 



a(p) 



(2.18) 
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Notice that this curve can be obtained from the generating function for disk amplitudes 
v(p), or equivalently from the D6 brane superpotential W(z). 

Finally, we recall that open topological string amplitudes on Calabi-Yau manifolds 
have an underlying integrality structure first obtained in [61] and further refined in [181156]. 
In the case of genus zero amplitudes with h holes, integrality leads to an expansion of the 
form 

oo 

F^zu ■ ■ ■ , z h ) = E dk ~*< w 4 Wl ■ ■ ■ 4 Wh Q d ^ (2-19) 

/3,Wi d=l 

where Hq w are integer numbers and give open string analogues of the Gopakumar-Vafa 
invariants. 

2.2 The annulus amplitude 

We start our discussion of open string amplitudes with the annulus amplitude, 

A(p,q)= E W*>~V, (2.20) 
k,e>i 

where k, I are winding numbers. It corresponds to the two-point planar amplitude in the 
matrix model, 

W (p,q) = d p d q A(p,q), (2.21) 

with the expansion, 

W (p, q)=Yl MFoM(t)P~ k ~ 1 q~ t ~ 1 - (2-22) 

k,t>l 

The amplitude Wo(p,q) is, together with the spectral curve itself, the most important 
ingredient in the analysis of loop equations, since it determines the rest of the amplitudes 
by recursion. In the case of the annulus amplitude it is also easy to make the connection 
between the theory of the chiral boson and the formalism of matrix models pETj. 

Let us then compute the annulus amplitude in the B-model, in the case of a local 
Calabi-Yau case described by (12. lip . As shown in [2], the B-model on (12. lip reduces to 
the theory of a chiral boson living on the curve E s associated to (12. lip . This chiral boson 
has a classical piece given by 

d(f) c i(u) =v(u). (2.23) 
The disk amplitude is given by the one-point function of this chiral boson, 

W{u) = (0(w)). (2.24) 
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Since only the classical part of the boson field contributes to this correlation function, one 
immediately recovers (I2.12p . The annulus amplitude should be given by the connected 
two-point function of the chiral boson, 

A(p,g)H<Kp)#7)> (c) - (2-25) 

Here, p, q are local coordinates on the Riemann surface chosen in such a way that, as 
p, q — > oo, A(p, q) vanishes. Only the quantum fluctuations of the chiral boson contribute 
to this connected correlator. It is well-known that this correlation function is a natural 
object on the Riemann surface £ 5 which can be expressed in terms of theta functions. 
Let us define the Abel map, 

rp 

Ui(p) = / Ui{p')dp, i = 1, • • • ,g, (2.26) 

•Jpo 

where Ui(p)dp is a basis of Abelian differentials on S 5 and p is a basepoint. We then 
have, 

A(p,q) =logtf*(u(p)-u(g)|r)-log(p-g). (2.27) 

where the subscript * means that the theta function has an odd characteristic, and we 
have subtracted the singular part of the correlator. The closely related quantity (I2.2ip is 
given by 

W {p, q) = mP)d<P{q)) ic) = B{p, q) - (2.28) 

where B(p,q) is sometimes called the Bergmann kernel. It can be also written as 

B{p,q) = d p d g logE(p,q), (2.29) 

where E(p, q) is the prime form on E 9 . Equation (12.271) generalizes the computation of 
the annulus amplitude for a target torus in [39], chapter 35.4, to mirror geometries given 
by a general Riemann surface £ 9 . The result ( I2.28P for the two-point connected function 
is a standard result in the matrix model literature, see for example [291 [32]. 

As it stands, the formula (12.271) is not easy to use in explicit computations. Fortunately, 
the Bergmann kernel B(p, q) plays an important role in the theory of matrix models since 
it gives an explicit expression for the correlation function (12.251) (see for example [291 E2] 
and references therein) . Its most important property is that it only depends on the spectral 
curve through the positions of the branch points X{. We will assume in the following that 
there are 2s branch points coming in pairs. The fact that Wo(p,q) only depends on Xi 
was proved in [13] for s — 1, in [5] for s = 2, and for general s in [11] . For s = 1 the 
spectral curve has genus zero and Wo(p, q) is given by [13] 

xxr, \ 1 ( pq- i(p + g)(si +x 2 ) + Xl x 2 \ 

Wo{p,q) = - ^ 7= = =, =, r " 1 )- (2- 30 ) 

2{p - q) 2 V y/(p - x x ){p - x 2 ){q - Xl ){q - x 2 ) ) 



S 



For s = 2 the spectral curve has genus one, and Wo(p,q) can be written in terms of 
Weierstrass function (see, for example, section 3.6 of [IS]). An explicit expression in 
terms of the branch points was found by Akemann in [8]: 




W (pq)- 1 ^ \{p-Xi){p-x A ){q-x 2 ){q-x 3 ) 



4(p - q) 2 V V (P~ x z)(P ~ x 3)(q ~ x\)(q - x 4 ) 



+ (P ~ x 2 )(p ~ x 3 )(q - Xi){g -x A ) \ , 231 s 
y (p - x^ip - x 4 )(q - x 2 )(q - x 3 ) _ 
(zi -x 3 )(x 2 -x 4 ) E(k) 1 



^UUp-^-^) K{k) 2{p - qr 

where K(k), E(k) are elliptic functions of the first and second kind with modulus 

k 2 = (Xl - X A )(X 2 - X 3 ) ^ 

(^1-^3)^2-^4)' 

It can be easily checked that these expressions have the expansion in p, q given in (I2.22p . 
One finds, for example, that 

F ,i,i(t) = (Xl 7/ 2)2 , F 0|1 , 2 (i) = l(x! - x 2 )\x 1 + x 2 ), (2.33) 
lo 62 

for s — 1. For s = 2 one has 

F ,i,i(t) = 7(^1 - x 3 )(x 2 - a; 4 )^77T + 77(^1 + a; 4 - x 2 - x 3 ) 2 . (2.34) 
4 K(k) lb 

In order to apply these formulae in the toric case, we write the spectral curve (I2.17f) from 
the disk amplitudes, and we read the position of the branch points from the zeros of the 
square root involved in the expression. How do we order them? The branch points come 
in pairs and define cuts in the complex plane. Each cut corresponds to a fixed "filling 
fraction." For example, for s = 2 the cuts are (xi,x 4 ) and (x 2 ,x 3 ) [25]. In the B model 
context, we will require that the cuts shrink to zero size at large radius/large complex 
structure, so that k 2 = at this point. This will lead to an expansion of the amplitudes 
compatible with the A-model interpretation at large radius. 

2.3 Holomorphic anomaly and Ray— Singer torsion 

The annulus amplitudes A(p, q), Wo(p, q) defined by (12. 27ft . (12.281) are holomorphic, but 
they are not modular invariant. Under a symplectic transformation 

t -> (At + B)(Ct + Dy\ (2.35) 
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the theta function transforms as 

i 

0„(u|r) -> (det(CV + 2 exp[i7ru(CV + D^Cu^uIr), (2.36) 

and Wo(p,q) transforms with a shift! 

W (p, q) -> W„(p, ?) " 27ri^(p)(Cr + D)- l Cuo{q). (2.37) 

Following the philosophy of [231 HQl EE], one can obtain a modular invariant object by 
adding an antiholomorphic piece in the closed string moduli, 

W (p,q) = W (p,q) - 7Tio(p)(lmT)~ l uj(q). (2.38) 

Notice that this object is still holomorphic as a function of the open string moduli p, q, 
and satisfies a very simple holomorphic anomaly equation for the closed string moduli. 
The Kahler metric on the Jacobian is 

9ki = ~K T - f )kh (2.39) 

and the tensor 

C ijk = Q -§ (2.40) 

is the three-point function in the closed string sector. The antiholomorphic dependence 
of W (p, q) is captured by 

dW {p,q) — y 

-= k = 2muJi{p)C- k ujjiq), (2.41) 

where we have raised the indices with the inverse metric g l K In the context of the 
holomorphic anomaly equation of [H] one can construct a propagator S 1 ^ satisfying 

d- k Sv = Cl, (2.42) 

therefore one can integrate (12.411) to write 

W (p, q) = 2mu i (p)S ij u; j (q) + f, (2.43) 

where / (the ambiguity appearing as an integration constant) is holomorphic in the closed 
string moduli. 



2 The modular properties of open and closed amplitudes in matrix models are extensively discussed in 

E2. 
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It was shown in [T4] that in the B-model the annulus amplitude should be related to 
the Ray-Singer torsion of the d operator coupled to a vector bundle V [66]. This bundle 
is associated to the D-branes that provide the boundary conditions [7T]. In fact, the 
Ray-Singer torsion should provide a natural modular invariant completion of the annulus 
amplitude (12.271) . We now discuss the relation between both quantities. The relevant 
background on differential operators on Riemann surfaces can be found in [TUl 124] . 

Flat line bundles over a Riemann surface S g are parametrized by the Jacobian J(H g ). 
Let us then consider the line bundle C associated to 

v = u(p) - u{q) E Jac(S 9 ). (2.44) 

In terms of divisors, this is the line bundle 

C=[p-q}. (2.45) 

We now consider the Cauchy-Riemann operator 

d c : n{£) ^Q '\C). (2.46) 

The Ray-Singer torsion of the 8 operator coupled to V = £ is simply the regularized 
determinant of the Laplacian Ac = d c dc- This has been computed for g = 1 in [66] by 
direct evaluation. For g > 1 one can regard the operator Be as a family parametrized 
by the Jacobian J(S S ), and the analytic torsion can be determined by using the Quillen 
anomaly. An explicit computation has been done in [42], but the result gets obscured by 
the fact that for g > 1 the dc operator has a nontrivial cokernel of dimension g — 1. In 
particular, the result of |42j depends on a choice of a point po G S ff . The Quillen anomaly 
gives a closed expression for 

\}p<M)? detA£ (247) 

det [XilXj) 

where Xi? i — 1> • • • , 9 ~ 1 is a basis for (Coker^)* ~ Ker d£-i® K , and W denotes a 
Wronskian. The appearance of the determinant of the inner products of the Xi is typical 
for computations which use the Quillen anomaly (see for example pi)])- The result of [42] 
for (12.471) is given by 



det'A 



Vol(E 9 )det Imr 



G(v,a,f3) = I u "7r 1 exp -27rImv(Imr) _1 Imv ■& " (v|r) (2.48 



(3. 



times a nontrivial function of po and r. In (12.481) . A is the Laplacian on scalars, and as 
usual the ' means that we remove the zero modes, (a, (3) is a characteristic induced by 
the choice of basepoint po through 

IC P0 = ra + f3, (2.49) 
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where JC Po is the vector of Riemann constants associated to the basepoint p$. The Quillen 
anomaly measures the failure of holomorphic factorization of the Ray-Singer torsion and 
reads 



d v d v 

V V det(xi\Xj) 



2mg dv{ A dv^. 



(2.50) 



For g = 1 the result for the Ray-Singer torsion is much simpler since the dependence 
on the basepoint p disappears, and one finds 



det A c = G(v, 1/2, 1/2) = exp -27rImv(Imr)" 1 Imv 



x) 



1/2 
1/2 



V T 



T}(t) 



(2.51) 



which indeed exhibits an odd characteristic. Therefore, for g = 1, the quantity 

A(p, q) = log det Ac — log \p — q\ 2 



(2.52) 



provides a modular invariant completion of (I2.27P which is antiholomorphic in both the 
closed moduli r and the open moduli v. The Quillen anomaly (12.501) gives a holomorphic 
anomaly equation for the open moduli. 

For g > 1, the relation between (12.271) and the Ray-Singer torsion is more complicated 
due to the presence of zero modes. If one chooses the basepoint p in such a way that 
the characteristic in (I2.48P is odd, then we can use the Ray-Singer torsion (or rather the 
quotient ( 12.471) ) as the natural completion of A(p, q). Notice that a natural generalization 
of (I2.52p for g > 1 is simply to set 



A(p, q) = log G(v, a, (3) - log \p - q\ 



(2.53) 



where G(v,a,/3) is defined in (I2.48P and we take (a,/3) to be an odd characteristic. This 
is a modular invariant quantity which satisfies the Quillen anomaly (12.501) and has (12.271) 
as its holomorphic piece. As shown in [TO], G(v, a, (3) can be regarded as the determinant 
of the Dirac operator on S 9 coupled to the flat bundle C, and for a choice of odd spin 
structure. This suggests that when considering D-branes on S 9 , the natural vector bundle 
V associated to them is not just C, but S^ 1 ® C, where S is the spin bundle associated 
to an odd characteristic. This is also consistent with the fact that D-branes behave as 
fermions \2\. 



2.4 Other topological string amplitudes from matrix models 

As we have seen, the annulus amplitude is particularly useful since it can be expressed 
simply in terms of the positions of the branch points of the spectral curve. The main result 
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of [T2j l8l [29l 120 ] 132] is that all free energies and correlation functions in the 1/N expansion 
of one and two-matrix models can be obtained recursively from the knowledge of the 
spectral curve and the two-point function (I2.28p . The recursion relation of [2"9"| I3"2" | |2"T | [52] 
involves evaluation of residues of meromorphic forms on the curve (II. ip at the branch 
points. In order to write down this recursion, we first notice that near any branchpoint 
%i there are two points q, q with the same coordinate x(q) = x(q). We now introduce the 
one-form, 



1 



9 



dE q (p) = - J d£B(£,p), (2.54) 

where B(£,q) is the Bergmann kernel (I2.29p . The recursion relation for the correlation 
functions reads, 



^ S dE (p) ( 
W g (p,p!,--- ,p h ) =y2~Res q=Xi -—— g,g ( W g -i(q, q,pi, ■ ■ ■ ,p h ) 

~ {y{q) - y{q))dx{q) \ 
9 h 1 \ 

+ Yl m \(h _ m l) W 9-e(<l>P«m> ■ ■ ■ ,Pa{m))W e (q, Pa{m+1) , ■ ■ - ,Pa(h))j 



(2.55) 



£=0 m=0 a&S{h) 



where the sum over a G S(h) is over permutations of h elements. There is also an 
explicit expression which computes the amplitudes F g in the closed sector from the open 
amplitudes W g (p). In that sense, open string amplitudes are more fundamental. The 
resulting formulae for the amplitudes require only the following information about the 
spectral curve: 

d k M(xi 
dp k 

i.e. the position of the branch points and the "moments" of the function M(p). In 
principle, the computation of these amplitudes is straightforward, but explicit expressions 
become cumbersome very quickly, specially when the number of cuts s > 1. We list some 
of them here, since they will be used in the tests performed in this paper. 

The three-point function at genus zero can be written for general s and was derived 
in [29]. Define first, 



x » 3 k ~ » * = ™, k>0, (2.56) 



G (p, q) = a{p)a{q) ( AW (p, q) + — — - ) . ( 2.57) 



Then, we have 



m / \ 1 G (x h p)G (x i ,q)G (x i ,r) 

W {p,q,r) = - > — . 2.5^ 

8a(p)a{q)a(r) M{Xi)a'(xi) 2 
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For s = 1 this simplifies to 

x 1 -x 2 / (p - - xi)(r - (p - x 2 )(g - x 2 )(r - x 2 ) 

^ ^ " 8(a(p)a(q)a(r)f 2 V M(^) M(x 2 ) 

(2.59) 

One can also write down results at higher genus. Knowledge of the first moments M(xi) 
makes possible to write down an explicit expression for F\. General formulae have been 
written down in [271 [T9] . For s = 1 we have the simple formula, 

F 1 = -^log[M(x 1 )M(x 2 )(x 1 -x 2 ) 4 }, (2.60) 

while for s = 2 we have [8j |4"5] 

1^1 1 
Fl = ~24 ln M fo') ~ 2 ln ~ 12 S ln ^ 4 ~ + 8 ln ^ 1 ~ X ^ + 8 ln ^ 2 ~ 

i=l i<j 

(2.61) 

One can also compute explicitly Wi(p), which involves the second set of moments M'(xj), 
in the cases s = 1,2, by using the results of [121 18], respectively. We emphasize that all 
these results are particular cases of the recursive procedure developed in [29], [201 E2] that 
we just reviewed. 

Although the above expressions were derived in the context of matrix models, we claim 
that they will lead to the correct open and closed B-model topological string amplitudes 
on a general local Calabi-Yau manifold, provided we use the spectral curve derived from 
the disk generating functional. In the next two sections we will give supporting evidence 
for our claim by studying some well-known toric CYs. 



3 Local curves 

3.1 Closed string sector and phase transitions 

By local curves we refer to the family of toric Calabi-Yau manifolds given by the total 
space of a bundle over a sphere, 

X p = 0{p-2)®0{-p) ^P 1 , peZ. (3.1) 

A notational warning: in this section, p will denote the integer labeling X p , and not a 
point on the spectral curve as in the rest of the sections of this paper. Topological string 
theory on X p has been very much studied recently (see [6] and references therein). As 
explained in [17], the A model has to be defined equivariantly. The most natural choice 
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(the equivariant Calabi-Yau case) corresponds to the antidiagonal action on the bundle. 
As usual, the topological closed string amplitudes F g (t) are put together into a single 
total free energy, 

oo 

Fx p (g s ,t) = J29 2 r 2 F g x ?(t) (3.2) 

3=0 

and its exponential is the closed partition function 

Z Xp = exp F Xp (g s ,t). (3.3) 

Here, t is the only complexified Kahler parameter of X p , and it corresponds to the size of 
P 1 . 

Closed and open partition functions on toric geometries can be computed by using 
the topological vertex introduced in [I], and we will rely extensively on it to test the 
B-model expressions. We collect here some formulae from the theory of the topological 
vertex which will be useful in the following. First of all, we define the q-number [n] as 

[n] = q n ' 2 - <T"/ 2 , q = e 9a . (3.4) 

A representation R of U(oo) is encoded by a Young tableau, labeled by the lengths of its 
rows {li}. The quantity 

l(R) = Y,U (3-5) 

i 

is the total number of boxes of the tableau. Another important quantity associated to a 
tableau is 

K R = ^2u(li-2i + l). (3.6) 



We also introduce the quantity 



oer 
as well as 



W R = q~^ J] 1 (3.7) 
- LJ - [hook(n)J 



W RlR2 = ? («« 1 +«h 2 )/ 2 W R{/R W R y R , (3.8) 



R 

>/ 



where R denotes the transpose tableau (i.e. the tableau where we exchange rows and 
columns), the "skew" W R is defined by 

w r/r/ = J2 n r'R" Wr " ( 3 - 9 ) 

R" 
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and Nft iR2 is the Littlewood-Richardson coefficient defined by the tensor product of U(oo), 

R 

Both (13.71) and (13.81) are particular cases of the topological vertex, which is given by 

Cr 1 r 2 r s — Q £ N q£ N q£ 3 2 W (3- 11 ) 

The topological string partition function on X p is given in terms of the Wr defined in 
Z x P = E W r Wr4 p ' 1)rr/2 Q 1{R \ Q = ("l) P e-*. (3.12) 

R 

Although (I3.12p gives an all-genus expression, it is effectively an expansion in powers 
of Q. In order to obtain an expression for each F g p (t) to all degrees in Q, one usually 
appeals to mirror symmetry. Since mirror symmetry in this equivariant setting is not so 
well-developed (see however [33] for recent progress), in [18j we analyzed the saddle of 
the sum over tableaux (13.121) directly. We found that the saddle is encoded in a density 
of tableaux p(A) with support on the interval 

(x 2 ,xi) U (xi,e), (3.13) 

where X\ i are nontrivial functions of the Kahler parameter t. To specify these, one 
introduces the mirror map 

The endpoints of the cut are given in terms of ( by 

Xl = (l-C)- p (l + C h ) 2 , x 2 = (l-C)- p (l-d) 2 - (3.15) 



We note that the variable ( is related to the variable w introduced in (18] by £ = 1 — w. 
This saddle-point solution was used in [18] to show that the above model has critical 
behavior for p > 2. At the critical point 

<-=^hp (3 - 16) 

the theory undergoes a phase transition in the universality class of pure 2d gravity. In 
fact, if one takes the double-scaling limit 

C -> Cc g s -»• 0, z fixed, (3.17) 
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where 



^ 5/2 = ^ 2 ^4^(Cc-C) 5 , (3.18) 



4(1 -Cc) 3 

then the total free energy ( 13. 2\\ becomes the free energy of 2d gravity 



F<m(z) = -U 12 ~ ^ log z + £ a^' 5(9 " 1)/2 , (3-19) 
where the coefficients a g can be obtained by solving the Painleve I equation 



u 2 - -u" = z (3.20) 
6 

satisfied by the specific heat 

u(z) = -F'^ ) (z). (3.21) 
3.2 Open string amplitudes 

We will now consider toric D-branes in the outer leg of the toric diagram associated to 
the geometry X p , as shown in Fig. [U We will focus here on p > 2. For p = 1 (the 
resolved conifold) there is an explicit matrix model description of the background (the 
Chern-Simons matrix model of [31]) and one can use this to study D-brane amplitudes 
[621 PTT] . It is well-known that open string amplitudes associated to non-compact branes 
need the specification of a framing [5j. Here, as well as in the rest of the examples, we 
will set the framing to zero. 

The open string amplitudes associated to toric D-branes can be computed in the A 
model by using the topological vertex, and they are encoded in the generating functional 

F{V) = log Z(V), (3.22) 

where 

Z(V) = —- J2z R Tr R V. (3.23) 

ZX * R 

In this example we have that 

Zr = J2 W RS q {p - 1)Ks/2 Q iiS) W s . (3.24) 
s 

The functional ( 13.231) is related to the generating functions (12. 4ft as follows: 



oo oo 



F(V) = J2J29 2 r 2+h F 9 (^--- ,z h ), (3.25) 

9=0 h=l 
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Figure 1: The outer brane in X p 



after identifying 



■ Z 



"'A 



TrV 



TrV 



(3.26) 



According to the general philosophy explained in the previous section, the first step 
in order to compute the open string amplitudes is the identification of the spectral curve. 
This can be obtained from the disk amplitudes. If we had a well-defined mirror symmetry 
setting for local curves we could derive the superpotential with the techniques of [7J [5]. 
However, here we do not know of any construction like (12.111) which allows to determine 
W(z) by integration, and we have to use another approach. Fortunately, the analysis of 
[T8] gives a detailed description of the saddle of the sum over partitions in terms of an 
effective matrix model, and we should expect that the resolvent of such a model is the 
object we need. Indeed, let us define 

p(v) 1 . A 



u (X) = / - 



e dv 
v 



X2 



7 log \ 
A A — e 



(3.27) 



\-v 

Notice that this differs from the resolvent as defined in [18] . but it is more natural since it 
integrates the density of tableaux over its full support. After some calculations one finds, 



a; (A) 



1 p 
A + 2A l0g 



a/(A - xi)(A - x 2 ) - A - y/xix 2 



+ A bg 



(vT 



Vx 



X2 



4A 



We can write this, using (I2.16p . in the form 

1 



u> (X) = -(V'(X)+y(X)), 



(3.28) 



(3.29) 



where 



V'(X) = ^\ogX- t - 



(3.30) 
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plays the role of the potential. The spectral curve is 



y(X) = - (tenli" 1 



y/{\- Zi)(A - X 2 ) 



A 



X1+X2 
2 



ptanh 



-1 



y/(X - Xl)(A - X 2 ) 

A + v /x 1 x 2 



(3.31) 



One can check that the resolvent (13.281) is the generating function for disk amplitudes. 
After using (12. 9p one finds, 



W{z) 



+ 



l + Q 
\ + Q 



l+p) pQ 2 + p (-l + 6p-8p 2 + 3p 3 ) Q 3 ^ 



2 6 
[l + 2p 2 )Q 2 p (-2 + 3p-4p 2 + 3p 3 ) Q 3 



4 



6 



+ ■ 



z 2 + - 



(3.32) 

This is precisely what one obtains in a topological vertex calculation (and leads to integer 
disk invariants for all p after taking into account the sign (— l) p in (13. 12ft ). This also 
confirms that (13.311) is the spectral curve for this model. 

As we explained in the previous section, all open and closed string amplitudes in this 
background can be derived from the knowledge of the branch points and of the moments 
of M(A), which involve the values of M(A) and its derivatives at these branchpoints. The 
branch points of the curve H3.31I) are of course x\,x 2 , and M(A) is obtained from the 
defining equation (II. ip . Using this information, together with the mirror map (I3.14p . we 
can compute the annulus amplitude from (12.301) . We obtain the expansion 



2 

A( Zi ,z 2 )=(q + P 2 Q 2 + ^ (l-4p + 3p 2 ) Q 3 
2 

+ — (1 -9p + 24p 2 -24p 3 + 8p A )Q 4 + ■■ 
3 

+ (q+{1+P + P 2 )Q 2 + ^ {-l + 2p + 3p 3 ) Q 



Z\Z 2 

3 



(3.33) 



+ ^ (-l + 8j9-20p 2 + 21j9 3 - 24p 4 + 16j9 5 ) Q 4 + ---)^ 2 2 + --- 

in agreement with a vertex computation for the outer brane. It is also easy to compute 
(12.601) . The first moments are given by 



M(x x ) = — 



P 



M{x 2 ) = — 



X 1 \X 1 -X 2 y/x~[(y/xl+ y/X2) 

One finds, up to an irrelevant numerical constant, 

'(P-1) 2 C(Cc-C) 



v 



X 2 \X 2 -X\ ^Jxl>(^/X^ + ^/X~2) J ' 

(3.34) 



Fx = log 

1 24 s 



(3.35) 
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This genus one amplitude can be written as 



Fi = F? at + log Q, 



(3.36) 



where F[ nst is the instanton part of Fi which was presented in eq. (6.1) of [18J and 
rederived in |33j. Of course, this reproduces the answer obtained from (13.121) . We have 
explicitly checked that the open string amplitudes with g = 0, h = 3 and g = h = 1 
derived in the matrix model formalism agree with vertex computations. Notice that the 
g = h = 1 amplitude involves the derivative of M(A) at the branch points, and therefore 
it probes in detail the structure of the spectral curve (13.311) . 

3.3 FZZT branes, instantons, and large order behavior 

As we reviewed above, for p > 2, topological string theory on X p leads to 2d gravity in a 
suitable double-scaling limit around the transition point in Kahler moduli space ( = ( c . 
We showed this in [18] by analyzing the closed string amplitudes. Now we have closed 
expressions for open string amplitudes (in particular, for the disk and the annulus), and 
we can analyze the double-scaling limit from the open string sector point of view. 

First of all, we recall that the critical point has a natural interpretation in terms of the 
spectral curve describing the model: it corresponds precisely to the value of the parameter 
C (or Q) for which a zero of the function M(A) coincides with the endpoint of the cut x\. 
This happens when 



which indeed holds at the critical point ( 13.161) . 

Let us now consider the double-scaling limit of the disk and the annulus amplitude. 
For the disk amplitude, the scaling part is precisely the spectral curve (I3.3ip . We first 
notice that near criticality the endpoints behave as 



M( Xl ) = 0, 



(3.37) 



Xl = 4 C) + e(c - c) + o(c - c) 2 , x 2 = 4 C) + o(c - c c ) 2 



(3.38) 



where 





9 = 2(1 



(1 



(1 



Q- 



p2 

(g - 2) 2 

{ P -ir 
P p(p - 1) 

p- 2 



(3.39) 




(3.40) 
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This defines the "renormalized" open string modulus s and says that, in the scaling limit, 
the brane must be located near the endpoint of the cut. We now reexpress (I3.3ip in terms 
of the scaling variables s and z, which is given in (13.181) . After some cancellations, we 
find that 

— y{X)d\ -> y(s)ds = -—z*(2s - l)VTT^ds, (3.41) 

9a 3 

for all p > 2. The r.h.s. is nothing but the Laplace transform of the macroscopic loop 
operator of 2d gravity [3U [25], and it corresponds to the disk amplitude of the FZZT 
brane of Liouville theory (see [671 [52] an d references therein). Notice in particular that, 
up to a normalization, the double-scaled spectral curve (I3.4ip is nothing but 

T 2 (y)=T 3 (s), (3.42) 

where T n (cos2;) = cos(nz) are Chebyshev polynomials of the first kind. This is the 
spectral curve of the (2,3) model [251 EI]- The main conclusion of this analysis is that 
the outer, toric D-brane of the local curve X p shown in Fig. [1] becomes, in the double- 
scaling limit, the FZZT brane of pure 2d gravity. One can further check this by taking 
the double-scaling limit of the annulus amplitude, which is 

W(X, //) - — ^ ( I , 2 + s + t = _ i) (3.43) 

up to an overall constant. Here we have expressed /i in terms of the renormarlized open 
string modulus t (not to be confused with the Kahler parameter). The r.h.s. of f !3.43|) is 
the Laplace transfom of the two-point function of the macroscopic loop operator [T3l 134"] . 
which corresponds to the annulus amplitude of the FZZT brane. Of course, the fact that 
the annulus amplitude H 2 . 3 1) has this limit when we scale its arguments in the way we 
have done it was noticed long ago in [13j . What is not obvious here is that annulus 
amplitudes for toric D-branes are given as well by an expression like ( 12.301) . This is a 
consequence of the formalism proposed in section 2, and it implies that the toric brane 
becomes a FZZT brane near ( = ( c . 

It can be easily shown from the Painleve I equation (13. 20f) that the coefficients a g of 
the perturbative genus expansion (13 . 19[) have the asymptotics 

a g ^A-^{2g)\ ) (3.44) 

where 

A = ^5, (3.45) 
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see for example [25jj. Based on general considerations in quantum field theory [49] we 
would expect that the large order behavior in (13.441) is controlled by instanton effects of 
the form 

exp(-AsS). (3.46) 

It was noticed in [681 [22] that this instanton has a natural interpretation in terms of 
eigenvalue tunneling. Recall that in the context of the matrix models, the variable x in 
(II. ip is an eigenvalue of the underlying matrix. The instanton effect is obtained as we 
move one single eigenvalue at the end of the cut x\ to a saddle point of the effective 
potential, located at the top of a potential barrier. This point xq satisfies 

M(x ) = 0, (3.47) 

and it is the zero of M(A) which collides with x\ at the critical point. The instanton 
action is given by 

A= / j/(A)dA. (3.48) 



XI 



As we approach the critical point, x\ — > xo and (13.481) vanishes. But in the double-scaling 
limit, since g s — * simultaneously, we find a nonzero result for A/g s . David explicitly 
showed in [22] that one obtains in this way the value (13.451) . Indeed, we can recover his 
result from our toric realization of 2d gravity. In terms of the renormalized open string 
variable s, the positions of xi, Xq are given by s = — 1 and s = 1/2, respectively, as one can 
immediately verify by looking at (13.401) and (13.411) . We can then compute the instanton 
action in the double-scaling limit to be 

- — [ 2 ds(2s- l)vTTs (3.49) 
3 J-i 

which agrees precisely with (I3.45p . 

We would like however to extend this analysis off criticality, in order to understand 
non-perturbative effects and the large order behavior of perturbation theory in topological 
string theory. In analogy with the 2d gravity limit, the relevant instanton effects will be 
obtained as we move the toric brane form the endpoint of the cut to a saddle point. The 
instanton action will now be a function of the couplings of the theory, as it is familiar 
from the analysis of instanton effects and large order behavior in the 1/N expansion (see 
[38] for examples in low dimensions). In our case, the action will depend on Q, and we 
will have a non-perturbative effect of the form 

exp(-^M), (3.50) 



3 The normalization in section 7.2 of this reference is appropriate for an even potential. The generic 
critical potential leads to an extra y/2 in A as compared to their computation. 
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leading to the large order behavior 

F 9 (Q) ~ (A(Q)r 2 °(2g)l (3.51) 

The instanton action A(Q) is given by the integral (I3.48P in the full theory, where y(X) is 
(13.311) and Xq is now the zero of M(A) which gives the dominant instanton contribution, 
i.e. the instanton with the smallest action in absolute value (recall that in general the 
instanton action might be complex). At least for Q near Q c , we expect that, due to the 
connection to 2d gravity, the point xo becomes the point s = 1/2 in the double-scaling 
limit. This means that it will behave like 

x = - \e{C - Q + 0(( - Cc) 2 (3.52) 

near the critical point. The solutions of M(A) = are easily obtained for any given p. 
For example, for p = 3, one finds a unique solution, 

4:X\X 2 

Xo = (jii-^r p=3 ' (3 - 53) 



which behaves indeed like (I3.52p . For p = 4 there are three solutions, and the one with 
the right critical behavior (I3.52j) is given by 



2-y /xTxo 

x = r- l= , P = 4 - 3.54 
\/X\ - Jx 2 



The integral ( 13.481) can be computed in closed form. One obtains, 

A(Q)=F(x Q )-F( Xl ), (3.55) 

where 

F( X ) = -log (/,(,)) (log (/,(,)) - 21og(l + { J M % r ) + >og(l + ( JH% )) 
OT . ( 2 fi(x) \ ( 2/i(x) \ (x 1 -x 2 ) 2 

- 2112 r ( ^ - V^) 2 > ~ 2 r ( v^r + v^W " log ~r- lQg - 



-plog(/ 2 (x))(log( /2 (,)) + 21og(l - - log(l - ^^0 ) 

- 2pLi 2 (--M^L) + 2pU 2 (—M^—) + ^(logx) 2 +plog(VxT+ v^) 2 log x 



2 y /x 1 x 2 J V( v /x7+ v /x 2 " 

(3.56) 

and 

/i0*0 = v 7 ^ -^i)(g -^2) +x - x 2 2 , 

f 2 {x) = y(x — Xi)(x - X 2 ) + X + y/XiX 2 . 



2.3 



Q>0 



Q<0 




X 2 



Figure 2: This figure shows the instanton effects in the p = 3 theory for real Q. When 
Q > the cut (xi,X2) is in the real axis and the instanton is obtained by moving one 
brane from x\ to xq. For Q < 0, x\$ are complex, with x\ = x\. The cut is parallel to 
the imaginary axis and there are two complex conjugate instantons obtained by moving 
the branes at the two edges of the cut. 

It is easy to check, for the values of x that scale like (I3.52p . that this action indeed 
becomes (13.451) in the double- scaling limit. 

Let us first analzye the instanton effects for Q real and positive in the range Q c > 
Q > 0. In this case, the branch points x±,X2, as well as Xq, are on the real axis, with 
< X2 < x\ < Xq, and the instanton effect corresponds to a single brane moving from X\ 
to Xq (see the left hand side of Fig. [2]). This leads to an instanton action which is real and 
positive and vanishes at the critical point Q c . In Fig. [3] we show the values of this action 
for p = 3,4 as a function of Q in this range. Since this instanton effect is of order l/g s , 
this analysis shows that, for p > 2 and Q real and positive, the genus expansion f)3.2p has 
the typical stringy divergence ~ (2g)! |68j and it is not Borel summable. 

It is interesting to compare this behavior to the one found in the Hermitian matrix 
model with a potential of the form 

^TrAf 2 + /tTrM 4 . (3.58) 

The F 9 (k) are analytic at k = with a radius of convergence r = 1/48. The genus 
expansion diverges as (2g)!, and the instanton action behaves very differently depending 
on the value of k. For — 1/48 < k < it is real and positive, and the genus expansion 
is not Borel summable. This of course reflects the instability of the model for negative k 
and leads to the non-Borel summability of 2d gravity. Topological string theory on X p 
with p > 2, Q c > Q > shows exactly the same behavior that the quartic matrix model 
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in the unstable region. 




0.02 0.04 0.06 0.08 0.1 



Figure 3: This figure shows the instanton action (13.551) as a function of Q c > Q > 0. The 
upper line corresponds to p — 3, while the lower line is p = 4. The action vanishes at 
Q c {p), and it is real and positive. 

Of course, in topological string theory the parameter Q is in general complex, due 
to the presence of the B field, and this will lead to complex instanton actions. Let us 
consider for example the case p = 3, and Q in the range — Q c < Q < 0, in which the 
genus expansion (HOD is still real and each F* 3 (Q) converges. For these values of Q the 
branch points x±,X2 are complex conjugate to each other, and the point xq moves to 
the negative real axis. The instanton effects are obtained by analytic continuation from 
the case Q > that we have just analyzed, where the instanton is real. The instanton 
obtained by moving one brane from x\ to xq is now complex, but since the free energy 
is still real, we should expect two instantons with complex conjugate actions. The other 
instanton is obtained by moving one brane from the other end of the Fermi surface, i.e. 
from X2, to x (see the right hand side of Fig. [2]). Since the instantons are now complex, 
the resulting genus expansion might be Borel summable. For general, complex Q, the 
cut (xi,X2) will be located in a generic direction in the complex plane, and one should 
analyze in detail the possible instantons obtained by moving the branes to the zeroes of 
M(A). The instanton actions in this more general setting can be computed as well from 
the analytic result (I3.56p . 

Our analysis indicates that, for real positive values of Q < Q c , the string perturbation 
series will bot be Borel summable, as in 2d gravity, but it also suggests that one can find 
Borel-summable expansions by rotating Q in the complex plane. As a final comment, we 
recall that instanton effects in 2d gravity can be interpreted in terms of ZZ branes [9]. The 
contribution of the ZZ brane can be written as the difference of two FZZT branes located 
at the points x\ and xq, as in ( 13.551) [671 E2]- The above computation in topological string 
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theory is very similar due to the underlying matrix model structure. 



4 Local surfaces 

In this section we discuss two of the most studied examples in topological string theory, 
namely the local Calabi-Yau manifolds obtained by considering the anticanonical bundle 
on P 1 x P 1 and on P 2 . We will make explicit the prescription of section 2 to compute open 
string amplitudes and present some examples. We will focus in all cases in outer branes 
with zero framing. 



4.1 Local 



x 



The local P 1 x P 1 geometry has two Kahler parameters Q s , Q t - The mirror can be 
described by the equation (12. lip with [5] 



F{e u , e v ) =e u + z t e~ u + e v + z s e~ v + 1. 



(4.1) 



The variables u, v are only defined up to a SL(2, Z) transformation, and one chooses one 
description or another depending on the brane that one wants to describe. The variables 
Zt,z s are the linear sigma model parameters, which are related to Qt, Q s through the 
mirror map. The explicit expression for this map is, 



Qt = z t exp 
Q s = z s exp 
The inversion of (14.21) gives 

z t = Qt ~~ 2QsQt ~ 2Qt ' 



2 E &L+*Z*lM 



k,l>0 



2 £ 

k,l>0 



(k\) 2 (l\) 2 

(2k + 21 - 1) 
(fc!) 2 (/!) 2 



(4.2) 



Q s - 2Q s Q t - 2Q\ + 



(4.3) 



Let us now consider the disk amplitude for an outer brane in this geometry, shown in 
Fig. HI This was computed in p2 ;50j5 The function v(z) is given by 



v(z) = log 



l-z + z s z 2 + 



z + z s z 2 ) 2 — Az t z 2 



(4.4) 



J This is brane III in the notation of 5\. The results presented in table 3 of this paper do not seem 



to be correct, however. A computation of the disk superpotential for the outer brane of local 
presented in example A. 2 of [50], and this can be checked against the topological vertex. 



is 
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Figure 4: The outer brane in local P 1 x P 1 . 

In order to make contact with enumerative results in the A model, we have to take into 
account the important fact that the open string modulus z is not a flat coordinate. There 
is a nontrivial open mirror map, and in the case of the outer brane in local P 1 x P 1 this 
is given by 

(4.5) 

where Z is the flat coordinate. The quantities that we compute in the B model will be 
functions of Zt, z s , z, and we will have to reexpress them in terms of the flat coordinates 
Qt,Qsi Z hi order to compare them with the quantities computed in the A model. 

As we explained in section 2, the spectral curve can be read from v(z). In terms of 
p = z~ l we find 

2 

y{p) = — tanh 
p 

This spectral curve has four branch points, given by 

xi,2 = \ - v ^T^V / ( 1 + 2 v / ^) 2 - 4 ^' *3,4 = - + v^±^(l-2v^) 2 -4^ s . (4.7) 

The labeling of the cuts has been chosen in such a way that the cuts {xi,x^), (#2,^3) 
shrink to zero size when Zt = z s = 0. We can now use the explicit expression ( 12.311) to 
extract the generating function A(p, q) of annulus amplitudes for outer branes, and we 
have to reexpress the result in terms of the flat coordinates. The open mirror map is the 
same for both the two open moduli, therefore we will set p = z^ , q = z^ 1 and use (14.51) 



p + p 



2\2 



4zp 2 



z s -p + p 2 



(4.6) 
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for both. One obtains the following expansion, 

A(Z U Z 2 ) = (Q t + 2Q t Q s + AQ 2 Q S + 6Q 3 Q S + 36Q 2 Q 2 S + 4Q t Q 2 s + 6Q t Q 3 s + • • • ) Z X Z 
+ (Qt + Qt + 1Q t Q s + 3Q t Q 2 s + 5Q t Q 3 s + 6Q 2 Q S + 30Q 2 Q 2 S + 10Q 3 Q S + ■ ■ ■ ) Z X Z 2 2 
+ (Qt + 1Q s Qt + ZQlQt + 4Q 3 a Qt + 3Q t 2 + 12Q s Q t 2 + A0Q 2 s Q 2 + Q 3 t + 2AQ S Q 3 + ■ ■ - ) Z 2 Z 
+ (Qt + 2Q s Q t + + 4 Q3Q t + M + 10Q S Q 2 + 33 q2q2 + Q 3 + 20Q S Q3 + . . .^2^ 

+ ••• 

(4.8) 

Notice that this satisfies the integrality conditions (12. 19[) . The above result can be com- 
pared to that obtained with the topological vertex. The total partition function for this 
geometry can be written as 

: £ Q? W(* 3 ) g « +W?Ei "^C^C^C^C^. (4.9) 



2 



2 v2 
2 



Ri 



M ( x i) = — Z. — i — 2V » = !,••-, 4. (4.11) 



The generating functional for the open string amplitudes in the presence of an outer brane 
is given again by (13.231) . where now is given by 

Z R = Y1 ctf^VQ***™**^ "^Cj^C^C^qC^. (4.10) 

Ri 

One can easily check that the annulus amplitudes computed from these vertex formulae 
agree with the expansion (14. 8p . 

It is possible to check as well the open, three-hole amplitude at genus zero and the 
closed, genus one amplitude F\, by using the formulae presented in section 2. We need 
the first moments of the function M(p), which are given by 

2 

One finds, for example, 
F (Z U Z 2 , Z 3 ) = (Q t + 2Q s Q t + 3Q 2 s Q t + 5Q 3 s Q t + 7Q 4 s Q t + Qt + QQsQt + 40Q 2 Q 2 

+ 189Q 3 S Q 2 + 10Q S Q 3 + 2\1Q 2 S Q\ + UQ s Qt + ■ ■ -)z x Z 2 Z z + ■■■ , 

(4.12) 

while for the genus one, closed string amplitude we find from ( 12.61}) the expansion 

F (0 Q- 1 hzO Qs Q * Q " Q ° Qt QsQt Q " Qt 2Q " Qt Q * 
UVs,W 12 6 12 18 24 6 3 2 3 12 

Q S Q 2 37 Q 2 Q 2 Q 3 2Q S Q 3 Q\ 
2 6 18 3 24 

(4.13) 
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The worldsheet instanton contributions in this quantity are indeed the correct ones, but 
the classical part — log Qt/12 is asymmetric in that it only includes the t parameter. This 
is probably due to the fact that the disk amplitude we start with is itself asymmetric in its 
treatment of Qt,Q s - In fact, it is nice that this asymmetry does not affect the instanton 
part of Fi, which is correctly reproduced. 

This ends our tests for local P 1 x P 1 . We will now briefly consider the case of local P 2 . 

4.2 Local P 2 

The local P 2 geometry has one Kahler parameter Q. The mirror can be described by 

F(e u , e v ) = -e v + e u + l- z t e 3u ~ v , (4.14) 

where Zt is the linear sigma model parameter. The mirror map is given by 

(3Jfe-l) 



Q = z t exp 



3£ 

L fc>0 



(k\y 



(4.15) 



and its inversion reads 

zt = Q + 6 Q 2 + 9 Q 3 + 56 Q 4 - 300 Q 5 + 3942 Q 6 + ■ ■ 



(4.16) 



Figure 5: The outer brane in local P 2 . 



Let us now consider the disk amplitude for an outer brane in this geometry, shown in 
Fig. This has been computed in [SI ED]- The function v(z) is given by 



v(z) = log 



1 + 2:+ a/(1 + z) 2 - 4z t z 3 



(4.17) 
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The open mirror map for the outer brane is given by 



(4.18) 



where Z is the flat coordinate. After setting z = p 1 , we find the spectral curve 



z i 
v = - tanh 
V 



p(p(p + !) 2 - 4^) 



P(P + 1) 

The branch points of the curve are x± — plus the roots of the cubic equation 

x(x + l) 2 - Az t = 0. 

In terms of 

£ = (1 + 54^ + 6^^ (1 + 27^)) 3 



they are given by 



3 3 



3 3 



X4 = 



(g ~ I) 2 
3£ ' 



(4.19) 

(4.20) 
(4.21) 

(4.22) 



where uj = exp(2i7r/3). 

We can now compute the annulus amplitude for the outer brane. We find the expan- 
sion, 

A(Z 1 ,Z 2 ) = (-Q + 4Q 2 -35Q 3 + 400Q 4 -5187Q 5 + ---)zT 1 zT 2 
+ (Q - 3 Q 2 + 24 Q 3 - 261 Q 4 + 3288 Q 5 + ■ ■ ■ ) Z X Z\ 
+ {-Q + 4 Q 2 - 25 Q 3 + 248 Q 4 - 2970 Q 5 + • • • ) Z X Z\ 



(4.23) 



+ (-Q + -^-20Q 3 + 194Q 4 - 2299Q 5 + ---)zT 2 zT 2 + --- , 

which again is in agreement with the results for the outer brane amplitudes obtained from 
the topological vertex. In this case, the total partition function is 



Zp2= (-I^^Q^'^Y^^C^C^Co^r,, 
and the amplitude for the outer brane with zero framing is obtained from 

Ri,R2,R 3 



(4.24) 



(4.25) 
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We can again compute the three-hole, genus zero amplitude. The first moments are given 
by 

M( Xi )= 2 z = l,...,4. (4.26) 

xfixi + 1) 

Notice that M(p) diverges at the first branch point X\ = 0. This is no problem for the 
three-hole amplitude f)2.58p . and simply means that the first term in the sum vanishes. 
One finds, 

F 3 (Z 1 , Z 2 , Z 3 ) = (-Q + 3 Q 2 - 36 Q 3 + 531 Q 4 - 8472 Q 5 + • • • ) Z 1 Z 2 Z 3 + • • • , (4.27) 

again in agreement with a vertex computation. Of course, for the expression f)2.6ip we 
should subtract the diverging moment for x% = 0. Once this is done, one obtains the right 
expression for F\. This ends the tests of the formalism proposed in this paper. 

A final comment on this example is in order. It is well-known that local P 2 has a 
conifold point where closed topological string amplitudes diverge. This point corresponds 
to the critical value 

4 C) = ~ (4.28) 

In terms of the spectral curve ()4.19|) . this is described by the collision of x 3 and x 4 , and 
it leads to a critical behavior that has been studied for example in [8} [3Q] - On the other 
hand, based on general properties of topological string theory, we should expect that the 
critical theory at this point is the c = 1 string at self-dual radius |36j. It was indeed 
suggested in [5HI EE] that the universality class associated to two colliding cuts should 
be indeed the c = 1 string, and some aspects of this correspondence have been further 
clarified in the series of papers [15]. It would be very interesting to analyze the open 
sector of topological strings on local P 2 near criticality, as we did in section 3 for local 
curves. This would make possible to extract nonperturbative effects due to D-branes and 
in particular to determine the large order behavior of the perturbative string expansion. 



5 Conclusions and open problems 

In this paper we have proposed a formalism inspired on matrix models which allows us to 
compute in closed form open and closed string amplitudes of the type B topological string 
on toric backgrounds. This formalism exploits the result of pfl, [2] which identifies the local 
B-model with the theory of a quantum chiral boson living on a Riemann surface. Since 
the general solution of matrix models on a spectral curve can be identified with such a 
quantum theory, we conclude that the rich results developed in the study of loop equations 
of matrix models can be applied in this context. 
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Notice that the "effective" matrix model description introduced here is rather different 
from the multi-matrix models of [53], 13]. Those models are parametrized directly by the 
flat coordinates, which appear either in the action or through the ranks of the matrix 
fields. In the formalism of this paper, the parameters entering into the description are 
not flat coordinates, and one has to use the mirror map (open and closed) to translate 
the answers into a suitable form for comparing with the A model. 

Of course, the results of Dijkgraaf and Vafa in [2H] provided a direct conection between 
matrix models and Calabi-Yau geometry, but the backgrounds analyzed in [28J are in a 
sense orthogonal to the ones that appear in Gromov-Witten theory and in the theory 
of the topological vertex. Here we have seen how to fit the latter in the matrix model 
formalism, and we have tested these ideas in various examples. The case of local curves 
is particularly interesting due to the connection of these backgrounds with 2d gravity. 
This makes possible to find the instanton effects responsible for the large order behavior 
of topological string theory on those manifolds. 

As a corollary of this matrix description we have a picture in which the toric Calabi- 
Yau geometry emerges as a large N description of a fundamentally discrete system of 
N interacting eigenvalues, and with a characteristic length of order 1/N ~ g s . This 
picture has already appeared in many systems which exhibit a gauge/gravity duality, in 
particular in the Dijkgraaf- Vafa scenario [28]. The results of this paper extend it to toric 
Calabi-Yau backgrounds. The discretization of the geometry that emerges is similar to 
that suggested in [61j . where a crystal model with characteristic length g s was proposed 
as the underlying stringy description of some toric geometries. It would be interesting to 
understand better the connection between these two pictures. 

Let us finish this paper with a list of open problems and directions for further research. 

• From a mathematical point of view, the results of this paper offer a wealth of mirror 
conjectures about the open string sector of toric Calabi-Yau manifolds. It should 
be possible to generalize the analysis of [37], which proved a mirror conjecture for 
the disk ampitude, to the case of annulus amplitudes and beyond. Of course, if 
our formalism holds in general, one should be able to prove direclty in Gromov- 
Witten theory that the generating functionals of open (or relative) Gromov-Witten 
invariants satisfy the recursion relations typical of loop equations found in [29] [32]. 
Since these invariants are given by sums over partitions, our formalism is equivalent 
to asserting that the corrections to the saddle-point of these sums satisfy the loop 
equations of an "effective" matrix model, at least for this class of examples. 

• More tests of our formalism are needed in order to understand amplitudes with 
higher g and h, and in particular the implications for the closed string sector. But 
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more importantly, in its current form this formalism only applies to amplitudes for 
outer branes, and one should generalize it to general configurations of branes in 
these geometries. A hint to understand this is given by the following observation: 
one can analyze the disk amplitude in terms of an open generalization of Picard- 
Fuchs equations [501 1511 135] . and it is known that outer and inner branes are on 
different sides of a phase transition which moreover involves different mirror maps. 
We expect that the same universal expressions in terms of non-flat coordinates (like 
(12.311) ) hold for all branes, but that one has to use different open mirror maps in 
different phases, leading in this way to the different expressions for open string 
amplitudes that one finds in the A model. After all, in the B-model outer and 
inner branes should be smoothly connected, since the Riemann surface described by 
(12.111) smooths out the toric diagram of the A modeM The mirror map for inner 
branes mixes open and closed coordinates in a nontrivial way [50], and one should 
probably work out a more general Picard-Fuchs system involving both the closed 
moduli and the h open moduli appearing in the general amplitude (12.41) . 

• As we mentioned in the introduction, it is not known what is the analogue of the 
holomorphic anomaly equations for open strings. However, by using the results 
of [32], it is possible [31] to find the holomorphic anomaly equations satisfied by 
the functions W g (zi, ■ ■ ■ ,z h ) defined by the recursive procedure of [291 [2H [32] (the 
holomorphic anomaly equation for the annulus amplitude is written down in ( 12.411) ). 
It is likely that these equations, which are in principle valid in the local case, give a 
strong hint of the structure of the holomorphic anomaly equations for open strings 
in the general case. Of course, one advantage of the procedure of this paper is that it 
is a recursive B-model calculation which does not exhibit the ambiguity of [H] . The 
holomorphic anomaly equations can be used to determine closed string amplitudes 
on local Calabi-Yau manifolds [46, 40], and it has been recently shown that for these 
backgrounds the holomorphic ambiguity can be fixed from a "gap" behavior at the 
conifold point [40] . It would be interesting to compare the formalism of this paper 
in the closed case with the procedure of [I6J HP] , and to explore in more detail the 
modular properties of the open string amplitudes. 

• From the results of [IE] and of section 3 of this paper, it seems that a way to 
understand certain nonperturbative aspects of topological string theory is to relate 
them to noncritical string theories, where these aspects are better understood. When 
discussing local P 2 we mentioned that there should be a double-scaling limit of the 
spectral curve (14.191) which gives the c = 1 string. A better understanding of 

5 This issue has been emphasized and discussed in [43] . 
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this limit would shed new light on nonperturbative aspects of topological string 
theory on local P 2 , like for example the identification of instanton configurations. 
In the case of P 1 x P 1 , there is besides the c = 1 limit a geometric engineering 
regime which leads to the SU{2) Seiberg-Witten theory [44J. The specialization of 
our formalism to this regime would give an effective matrix model description of 
Nekrasov's instanton expansion and its gravitational corrections [60]. That such a 
matrix model description should exist has been also pointed out in [i~5] . 
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